Energy avalanches in a rice-pile model 
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We investigate a one-dimensional rice-pile model. We show that the dis- 
tribution of dissipated potential energy decays as a power law with an 
exponent a = 1.53. The system thus provides a one-dimensional example 
of self-organized criticality. Different driving conditions are examined in 
order to allow for comparisons with experiments. 
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1 Introduction 



Statistical-mechanical investigations of driven nonequilibrium systems is a 
field of much current interest. A special class of such systems are those which 
reach the steady state through a self-organizing dynamics consisting of ava- 
lanches propagating through the system. If, in addition, the steady state is 
characterized by a power-law distribution for the sizes of avalanches the be- 
havior is referred to as self-organized critical (SOC) [l|]. Simple "sandpile" 
models have been introduced in order to illustrate the SOC behavior flrf|. 
Real sandpiles, however, have turned out to display noncritical behavior — i.e., 
they are described by avalanche size distributions with a characteristic scale 

Recently, an experiment with rice grains was performed |15| . It was found that 



a pile of elongated rice grains evolved into a SOC state with the distribution 
of dissipated potential energy E being flat for small avalanches and crossing 
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over to a power law with exponent a cxp ~ 2.0 for large avalanches. Physically 
motivated rice-pile models have been introduced which indeed do display SOC 
behavior in one dimension |I6HT9| 
that observed experimentally. 



but with an exponent which differs from 



In this paper, we measure the energy avalanches for the rice-pile model which 
we introduced in Ref. fl6|]. Previously we found that the distribution of av- 
alanches defined as the number of topplings s followed the power-law form 
P(s,L) ~ s~ T f s (s/L u °), with r = 1.53 and v, = 2.20 ||. Here, the distribu- 
tion P(E, L) of energy avalanches in a system of size L is measured. We find 
that P(E, L) scales with the same exponent as P(s, L) in accordance with 
the fact that a toppling event on average dissipates a constant energy. We 
change the drive in our model in order to make comparisons with the exper- 
imental results in \TE\ more quantitative. Despite the fact that this leads to 
a flat distribution for small E, in perfect agreement with the experiment, the 
power-law tail for large E is still described by an exponent 1.53. 



2 Rice-pile model 



We briefly describe the model we introduced in [[Lql . We consider a one- 
dimensional system of length L with a wall at i — and an open boundary at 
i — L + 1, where particles fall off the pile. The dynamics of the model consists 
of deposition and relaxation: At each time step one grain is added at i — 1. 
Then, the pile is allowed to relax in order to reach a new stable configuration. 
The relaxation process is considered to be fast compared to the deposition 
time scale. During the relaxation active columns topple one grain from % to 
i + 1 with probability p(5hi), where Shi = h(i) — h{i + 1) is the local slope. 
A column i of the pile is said to be active if, in the anterior time step, it (i) 
received a grain from column % — 1, (ii) toppled a grain to column i + 1, or (iii) 
column i + 1 toppled one grain to its right neighbor. The probability p(5hi) 
to move a grain is taken to be: 



p(Shi 



0, Shi < Si, 

p, Si < Shi < S 2 , 

1, ^2 < Sh h 



where we will use Si = 1, 5*2 = 4, and p = 0.6; the results are insensitive to 
the precise values of these parameters. 

Physically the parameter p describes the friction between rice grains. It also 
incorporates the possibility that a metastable packing configuration will be 
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reached when a grain topples. The friction effect is the new ingredient in the 
model compared to other models and it introduces a large range of slopes in 
the rice pile instead of a single critical value. The parameter Si accounts for the 
fact that small slopes are stable. The parameter S2 models the effect of gravity 
on the packing configurations. We assume that above the maximum value S2 
of the local slope, it is no longer possible for a local stable configuration to 
be achieved, thus a grain must be toppled. In the limiting cases p — 0, 1, or 
S2 = Si, we recover the model in Ref. [|l| (which has trivial behavior in one 
dimension) . 



3 Avalanche results 



We study the model in the slowly driven limit where the rate of deposition is 
slow enough that any avalanche, that might be started by a deposited grain, 
will have ended before a new grain is deposited. The simulations of the model 
show the existence of a SOC steady state. We follow the definition of Ref. 
15 1 and calculate the size of avalanches as the dissipated potential energy in 



between snapshots of the profile. Figure |I|(a) shows the probability density of 
avalanche sizes for different system sizes. We find that the scaling form 

P(E,L)~E- a f{E/L") (2) 



describes the distribution of avalanche sizes. The validity of Eq. (fj) is reassured 
by the good data collapse displayed in Fig. |l](b), where we used the values 
a = 1.53 ± 0.05 and v = 2.20 ± 0.05. Since a toppling on average dissipates a 
constant energy, we expect that a = r (and v = v s ) as observed numerically. 
By assuming that the average value is (E) ~ L in the critical state, it follows 
from Eq. (fj) that a = 2 — 1/V |2(J . This relation is in nice agreement with our 



numerical values for the exponents. We note that the scaling function / has a 
peak for values of the argument close to the cutoff region. This is a finite-size 
effect which is due to the possibility to form a supercritical state which then 



relaxes through a very large avalanche [16 



Our estimate for the exponent a is different from the value a exp ~ 2 reported 
in jy|, and the scaling function is also quite different from that observed 
experimentally. Part of the reason for this disagreement is due to the fact 
that the experimental conditions were different from how we measure the 
avalanche distribution P(E, L). The next step is thus to study the model under 
conditions as similar as possible to the experimental ones. In the experiment 
grains were deposited randomly in time with a rate of 4-6 grains between 



successive snapshots of the pile taken every 15 seconds [15]. Furthermore, since 



the profile of the pile was not known at every time instant when grains were 
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Fig. 1. (a) Log-log plot of the probability density P(E,L) for different values of L. 
The results was obtained with the parameters p = 0.6, Si = 1, and S% = 4. The 
data follows a power law distribution for several decades with a cutoff that depends 
on the system size, (b) Data collapse of the curves displayed in (a) according to 
Eq. (||) with the exponents a ~ 1.53 and v ~ 2.20. 
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deposited, the potential energy of the deposited (4-6) grains was estimated 
from the last profile obtained. The uncertainty in the added potential energy 
induces a noise level proportional to the system size (see below). 

To study the model under the "experimental conditions" we made the follow- 
ing assumptions: (i) snapshots of the profile of the pile are taken every N s 
time steps, and (ii) new grains are deposited on the pile at an average rate of 
1/Nd- Thus, every time step there is a probability 1/Nd of a new grain being 
deposited. In our simulations, characteristic values for iV s and Nj, were 10 4 
and 2100, respectively. 

Figure Q(a) shows the probability density of avalanche sizes for an accurate 
calculation of E. No flat part is observed for small values of E and no de- 
pendence on L is detected — except for finite-size effects. However, when we 
estimate E as it was done in the experiment of Ref. |HJ a significant change 
occurs. As shown is Fig. 0(b), a plateau whose height depends on L, and origi- 
nates from the uncertainty in the added potential energy, is observed for small 
E. For large E, there is a crossover to a power-law behavior at a value that 
scales with L. In fact, the data resemble quite well the experimental results. 
We find that the data in Fig. 0(b) are well described by the scaling form 

P(E,L)~L-P g(E/L»), (3) 

where g is a scaling function which is constant for small arguments and decays 
as a power law for large values of the argument. As can be seen in Fig. 0(c), 
a good data collapse is obtained with the exponents = = 1.00 ± 0.05 
for L < 640. On the other hand, for the two larger system sizes we again 
observe data collapse but to a shifted curve. The reason for this has to do 
with the possibility of an avalanche running beyond the time at which profiles 
are measured or a new grain is added: for small L almost every avalanche 
runs its course before a new grain is deposited or the profile is measured, but 
for large L we can still have running avalanches before a new picture of the 
profile is taken implying that large avalanches are not sampled and the entire 
distribution is shifted upwards for large energies. This hypothesis is confirmed 
by the fact that for the largest values of L the size of the biggest avalanche 
does not seem to grow. 

Another observation is that for the inaccurate estimation of the energy, the 
value of the exponent seems to change from 1.53 to about 1.4 [Fig. @(b)]. 
The reason for this can be understood as a simple finite-size effect and that 
it is only for a small region where the true power-law behavior is observed. 
One way to circumvent the finite-size effect is to measure instead the slope of 
the envelope which is close to 1.53. Thus, we expect the correct value of the 
exponent a to be 1.53. 
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4 Conclusions 



We investigate energy avalanches in a one-dimensional rice-pile model. We 
find that Eq. (|2J) provides a good description of our numerical results for the 
model when we drive it slowly. Furthermore, when we drive the model in a 
way close to the experimental conditions the numerical results are described 
by Eq. (|D in nice agreement with the experimental results in ||15|| . However, 
the numerical value for the exponent describing the power law decay is 1.53, 



6 




10° 



E/L 



Fig. 2. (a) Log-log plot of the probability density of avalanche sizes for different 
systems sizes. The data was obtained for conditions similar to the experimental 
ones but with an accurate determination of the dissipated energy. A power law 
dependence is observed with an exponent around 1.5. (b) Data from the same run 
as in (a) but now calculating the dissipated energy according to the method used 
in the experiment. It is remarkable how different the data looks from (a). A plateau 
for small E is visible and a dependence on L is detectable, (c) Data collapse of the 
curves shown in (b) according to Eq. (0) with the exponents = p, ~ 1. 



whereas the experiment gave the value a cxp ~ 2. Despite the fact that the 
overall scaling form can be understood by the rice model studied here, the 
disagreement in the values for a shows that the model needs extensions. One 
such possibility would be to include the effect of the kinetic energy of the 
particles in addition to the friction parameter p considered here. 
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